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Abstract. We present evidence for an expanding super- 
bubble in M82 (diameter: «130pc, expansion velocity: 
«45 kms" 1 , mass: rj 8 • 1O 6 M ). It is seen in the 12 CO 
(J = 1 — > 0), 12 CO(J = 2 1), 13 C0(J = 1 -> 0) 
and C 18 0( J = 1 — ► 0) lines. The superbubble is centred 
around the most powerful supernova remnant, 41.9+58. 
The CO observations show that the molecular superbub- 
ble already broke out of the disk. This scenario is sup- 
ported by ROSAT HRI observations which suggest that 
hot coronal gas originating from inside the shell is the 
main contributor to the diffuse X-ray outflow in M 82. We 
briefly discuss observations of the same region at other 
wavelengths (radio continuum, optical, HI, X-rays, ion- 
ized gas). From our spectral line observations, we derive a 
kinematic age of about 10 6 years for the superbubble. Us- 
ing simple theoretical models, the total energy needed for 
the creation of this superbubble is of order 2 x 10 54 ergs. 
The required energy input rate (0.001 SNyr -1 ) is reason- 
able given the high supernova (SN) rate of « 0.1 SNyr -1 
in the central part of M 82. As much as 10% of the energy 
needed to create the superbubble is still present in form of 
the kinetic energy of the expanding molecular shell. Of or- 
der 10% is conserved in the hot X-ray emitting gas emerg- 
ing from the superbubble into the halo of M 82. This newly 
detected expanding molecular superbubble is believed to 
be powered by the same objects that also lie at the origin 
of the prominent X-ray outflow in M82. It can therefore 
be used as an alternative tool to investigate the physical 
properties of these sources. 
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1. Introduction 

M82 is the best studied nearby starburst galaxy (D = 
3.25 Mpc). The central few hundred parsecs of this galaxy 
are heavily obscured by dust and gas which hides the cen- 
tral starburst region against direct observations at optical 
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wavelengths. Evidence for strong star-forming activity in 
the central region comes from radio (e.g. Kronberg et al. 
[1981 ) and infrared observations (e.g. Telesco etal. 1991) 
and also from the prominent bipolar outflow visible in Ha 
(e.g. Bland & Tully [1988| , McKeith et al. |1995| , Shopbell 
& Bland-Hawthorn 1998) and in X-rays (e.g. Bregman 
etal. 1995). The massive star formation (SF) is believed 



to be fuelled by the large amount of molecular gas which 
is present in the centre of M82. 

On the other hand, SF effects the distribution and 
kinematics of the surrounding interstellar medium (ISM). 
Recent mi llimet re continuum observations (Carlstrom & 
Kronberg 1990 ) suggested that the HII regions in M82 
have swept up most of the surrounding neutral gas and 
dust into dense shells. This is in agreement with the stan- 
dard picture: shells are created by young star-forming re- 
gions through strong stellar winds of the most massive 
stars in a cluster and through subsequent type-II super- 
novae (e.g. Tenorio-Tagle & Bodenheimer 1991 ). These 
processes are thought to blow huge cavities filled with 
coron al gas into th eir am bient ISM (e.g. Cox & Smith 
1974, Weaver etal. 1977). This hot interior is then be- 



lieved to drive the expansion of the outer shell of swept-up 
material. 

Once superbubbles reach sizes that are comparable to 
the thickness of a galaxy's disk, the bubble will eventu- 
ally break out into the halo. This then leads to an outflow 
of the hot gas with velocities much higher than the ex- 
pansion of the shell within the disk of the galaxy. In the 
following we present evidence for a molecular superbub- 
ble in M 82 which already broke out of the disk and seems 
to contribute significantly to the well-known prominent 
outflow of M82. 

2. Observations 

2.1. Molecular Lines 

For our analysis we used the 12 CO(J = 1 — ► 0) data 
cube obtained by Shen & Lo ( p.995[ ) with the BIMA ar- 
ray (spatial resolution: 2.5") an d the 13 CO(J = 1 — >• 0) 
data cube from Neininger etal. ( [1998 ) observed with the 
Platau de Bure interferometer (PdBI) (spatial resolution: 
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4.2"). In addition, we used unpublished PdBI data of the 
12 CO(J = 2^1) and C 18 O(J=1^0) transitions. These 
observations where carried out in April 1997 in the CD 
configuration, resulting in a spatial resolution of 1.4" x 1.2" 
( 12 CO(J = 2^1)) and 3.7" x 3.5" (C 18 0( J= 1 ->0)), and 
a velocity resolution of 3.3 kms" 1 and 6.8 kms" 1 , re- 
spectively. In order to increase the sensitivity to extended 
CO emission we combined the 12 CO(J = 2^1) data cube 
with single dish measurements obtained with the IRAM 
30m telescope. These observations were carried out in May 
1998. The combination is essential for this particular study 
because the receding part of the superbubble is only mar- 
ginally visible in the mere interferometer maps. A full ac- 
count of the data reduction will be given elsewhere. 



2.2. Evidence for an Expanding Superbubble 

Fig. shows th e inte grated 12 CO line emission published 
by Shen & Lo (1995). The cross corresponds to the po- 
sition of the supernova remnant 41.9+58 (SNR 41.9+58) 
which is the strongest cm continuum point-source in M 82 
(Kronberg et al. 1981). It is considered to be the aftermath 
of a 'hypernova', which exhibits a radio luminosity 50-100 
times greater than typical for type-II SNe (Wilkinson & de 
Bruyn 1990p . The line along the major axis indicates the 
orientation of the position-velocity (pv) diagrams shown 
in Fig. |. 
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Fig. 1. Integrated 12 CO line emission from Shen & Lo 
(1995). The lines along the major and minor axis of 
M 82 indicate the orientations of the pv-diagrams shown 
in Figs. | and | The cross marks the position of SNR 
41.9+58. 



The pv-cut is orientated in such a way that the signa- 
ture of the expanding superbubble is visible most dis- 
tinctly. The angular axis of the pv diagrams correspond 
to the offset in arcseconds from SNR 41.9+58. Besides a 



constant velocity gradient (cf. Shen & Lo 1995) the pv- 
diagrams show an expanding ring-like feature centred on 
SNR 41.9+58, with a central velocity of Vi sr « 150 kms -1 . 
The approaching velocity component of the ring is clearly 
seen in all cubes and is centred at 100 kms -1 . The reced- 
ing component is only marginally visible in the 13 CO 
(J = 1 -► 0) and C 18 0(J = 1 ->• 0) lines; its central ve- 
locity is about 190 kms -1 . An enlargement of the pv- 
diagram along the major axis is shown in Fig. |^ (left). 



From this diagram we estimate the radius of the ring to be 
(65 ± 5) pc and the expansion velocity (45 ± 5) kms -1 . A 
pv-diagram along the minor axis centred on SNR 41.9+58 
is presented in Fig. ^ (right). The orientation of the cut is 
shown in Fig. |l|. 



100 - 




offset 



Fig. 2. The pv-diagrams along the major axis of M82 
(orientation shown in Fig. [l]) for the 12 CO(J = 1 — ► 0), 
12 CO(J = 2 1), 13 CO(J = 1 0) and C 18 0(J = 1 ->0) 
cubes. The slice is centred on SNR 41.9+58. 



Fig. |^ (right) reveals two emission features centred at the 
position of SNR 41.9+58, which correspond to the ap- 
proaching (v rs 100 kms -1 ) and the receding component 
(y w 190 kms -1 ) of the expanding superbubble. Hardly 
any CO emission with a velocity between v= 100 kms -1 
and v = 190 kms -1 is found south and north of SNR 
41.9+58. The pv-diagrams therefore show that the ex- 
panding molecular shell has already broken out of the disk 
and now only shows the signature of an expanding molec- 
ular ring. It should be noted that the remaining molecular 
gas in M 82 shows clear solid-body rotation. 

2.3. Other wavelengths 

Radio continuum observations: 

Radio continuum observations at 408 MHz by Wills et al. 



(1997) unveiled a prominent 'hole' of approximately 100 
pc diameter around SNR 41.9+58 which they attribute 
to free-free absorption. They propose that this feature is 
due to absorption by a large HII region that has been 
photoionized by a cluster of early-type stars of which the 
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progenitor of SNR 41.9+58 was originally a member. It 
should be noted, however, that SNR 41.9+58 is only about 



50 years old (Wilkinson & de Bruyn 1990) and therefore 
cannot be the source that drives the expansion of the su- 
perbubble. But its presence supports the scenario that the 
expanding molecular superbubble is powered by an inte- 
rior stellar cluster. 
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X-ray observations: 

Fig. |^ shows an overlay of archival ROSAT HRI X-ray 
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Fig. 3. The pv-diagrams for the 12 CO(J= 1 — >0) transi- 
tion along the major and minor axis of M82. The spatial 
axis for each diagram is centred on SNR 41.9+58. The pv- 
diagram along the major axis is an enlargement of Fig. ^| 
(top). The orientation of the cut along the minor axis is 
indicated in Fig. [l]. 

Neutral hydrogen (HI): 

H I emission line studies show that the central kpc of M 82 
is dominated by absorption (Yun et al. 1993| ), which makes 
it difficult to study the H I kinematics in the central part 
of the galaxy. Recent H I absorption st udies against super- 
nova remnants in M 82 by Wills et al. ( 1998 ) disclosed two 
absorption features against SNR 41.9+58. The velocities 
of the components are 87 + 7 kms -1 and 200 ± 7 kms -1 . 
The HI component at 87 kms -1 can be attributed to 
the approaching component of the expanding superbub- 
ble. The absorption feature at 200 kms -1 , if associated 
with the receding part, seems to contradict the hypoth- 
esis that SNR 41.9+58 is located within the expanding 
superbubble. However, the velocity resolution of the ab- 
sorption study was rather poor (26.4 kms -1 ) and other 
explanations cannot be ruled out at this point. 



data (Bregman etal. 1995) onto the CO emission. To 



emphasize the diffuse, extended X-ray outflow we sub- 
tracted a model of the brightest X-ray point source at 
a = 9 h 55 m 50 s .4, 8 = 69°40'47.5" (J2000.0) and smoothed 
the residual image to a 6" resolution. The positions of the 
three point sources in the field are marked with stars, SNR 
41.9+58 is marked by a cross (see also Bregman etal.). 
The overlay shows that most of the diffuse X-ray emission 
arises from the vicinity of SNR 41.9+58. The shape of the 
diffuse emission is elongated, extending several arcminutes 
along the minor axis of the galaxy. Model calculations by 
Bregman et al. show that the emission is consistent with 
an outflow of heated material from the central region co- 
incident with the position of the expanding molecular su- 
perbubble. 
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Fig. 4. Integrated 12 CO(J = 1 ->0) line emission of M82 
(greyscale). The contour map represents the diffuse X-ray 
emission without the strongest point-source in the core of 
M 82 as observed by the ROSAT HRI. The three X-ray 
point -sources are marked with a star, the position of SNR 
41.9+58 by a cross; the centre of M82 (2.2 /zm peak) is 
indicated by a filled box. 



Optical observations: 

Optical observations of the centre of M 82 suffer heavily 
from light absorption by the prominent dust lanes, ren- 
dering an optical analysis of this particular region impos- 
sible. HST V- and I-band images of the centre of M 82 



(O'Connell etal. 1995) do not reveal any optical sources 
in the vicinity of SNR 41.9+58. The prominent Ha out- 
flow is visible north and south of the absorbing dust lane 
(Bland & Tully [l988| McKeith et al. ^995] Shopbell & 



Bland-Hawthorn 1998). The orientation of its filaments 



suggests that the Ha outflow emerges at least partly from 
the location of the molecular superbubble. 



We estimate an upper limit for the energy still present in 
the hot gas emerging from the superbubble. For the calcu- 
lation we adopt a density n = 0.2 cm -3 and a temperature 
kT = IkeV for the central region (Bregman etal.) as an 
upper limit for the mean values in the outflow. We find an 
upper limit for the energy of E gas < 2.7 x 10 53 ergs present 
in a cylindrical volume with r = 130 pc and z = 550 pc 
centred on the superbubble. From the thickness of the 
molecular disk of M 82 (w 90 pc) and the expansion veloc- 
ity of the superbubble we estimate a time elapsed since 
outbreak of r w 3 x 10 5 yr. The required outflow veloc- 
ity to reach the observed height above the disk therefore is 
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v » 1000 km s . This value is comparable with estimates 
based on particle aging as derived from radio continuum 



observations (Seaquist etal. 1985) and measured outflow 
velocities in Ha (McKeith etal. 1995). 



Ionized gas: 

Most studies of the ionized gas in M 82 unfortunately do 
not provide sufficient spatial and spectral resolution to re- 
veal details of the region under study. A high-resolution 
study of the Nell emission line has been published by 
Achtermann & Lacy ( 1995] ). The integrated line emis- 
sion in the south-western region of M 82 shows a similar 
double-peaked feature centred on SNR 41.9+58 as traced 
by the CO. The pv-cut along the major axis clearly shows 
a disturbed velocity field in the vicinity of SNR 41.9+58. 
The approaching side of the expanding superbubble at 95 
kms -1 is clearly visible while the receding side at 190 
kms -1 is only marginally seen. The diameter of the Nell 
shell (wlOOpc) seems to be somewhat smaller than the 
corresponding CO feature, suggesting that the inner part 
of the expanding superbubble is ionized. Similar features 
are visible in the distribution and the kinematics of the 



H41a emission line (Seaquist et al. 1996 ). The Ne II distri- 
bution might represent the transition from the hot coronal 
gas with a temperature of several 10 6 K towards the cold 
molecular rim of the superbubble. 

3. Discussion 

The picture that emerges from the observations presented 
above is the following: a major SF event at the centre of 
what today shows up as the prominent expanding molec- 
ular superbubble created a cavity filled with coronal gas 
by the combined effects of strong stellar winds and SN 
explosions. The pressure of the hot-gas interior drove the 
expansion of the shell of swcpt-up material until it broke 
out of the disk. The subsequent outflow of hot material 
today shows up as diffuse X-ray emission and contributes 
to the prominent Ha filaments. 

From the observational parameters of the expanding 
superbubble (expansion velocity and diameter) we esti- 
mate a kinematic age of the superbubble of 1 x 10 6 years. 
The kinematic age is an upper limit for the actual age 
since the superbubble is most probably decelerating. We 



use Chevalier's equation (Chevalier 1974) to derive the 
amount of energy needed to create the expanding super- 
bubble. Since Chevalier's equation applies to HI shells 
only we corrected the energy input by a factor of 2 to 
correct for the difference in mass between H2 and H. We 
estimate an ambient H2 density prior to the creation of 
the shell of order 120 cm -3 . This is done by converting 
the 12 CO(J = 1 — ► 0) line integral to H2 column density 
using a conversion ratio of N( H 2 )/W (CO) = 1.2 x 10 20 
cm~ 2 K _1 km _1 s (Smith etal. 1991) ) and estimating the 
volume of the region from which material was swept up to 
be cylindrical with z = 100 pc and r — 65 pc. We derive a 



total energy of order 2 x 10 ergs, which corresponds to 
an energy equivalent of approximately 1000 type-II SNe 
and the strong stellar winds of their progenitors. With the 
estimate for the age of the superbubble this leads to a SN 
rate of 0.001 SN yr _1 for the central stellar cluster. This is 
a reasonable number given the fact that the SN rate in the 
central part of M 82 was estimated to be about 0.1 SN yr" 1 



(Kronberg et al. 1981). The total H 2 mass of the superbub- 
ble is about 8 x 10 B M©. Less than 15% of the total energy 
is still present in the hot gas which emerges from the su- 
perbubble; the fraction of kinetic energy of the molecular 
superbubble is about 10% (E kin w 1.6 x 10 53 ergs). It 
should be noted, however, that the numbers given above 
are only order of magnitude estimates. 

The wealth of observations presented above indicates 
that the region under study represents an unusually ac- 
tive segment of the starburst in M82. The finding of an 
expanding molecular superbubble in this particular region 
supports this view and provides clear evidence for a foot- 
print of violent SF in the ISM of M82. Furthermore the 
analysis of the molecular superbubble presents an alter- 
native tool to investigate the energy release of the central 
source which drives the X-ray and contributes to the Ha 
outflow. 
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